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Abstract - A new, deadbeat type of direct torque control is 
proposed, analyzed and experimentally verified in this paper. 
The control is based on stator and rotor flux as state variables. 
This choice of state variables allows a graphical representation 
which is transparent and insightful. The graphical solution 
shows the effects of realistic considerations such as voltage and 
current limits. A position and speed sensorless implementation 
of the control, based on the self-sensing signal injection 
technique, is also demonstrated experimentally for low speed 
operation. The paper first develops the new, deadbeat DTC 
methodology and graphical representation of the new algorithm. 
It then evaluates feasibility via simulation and experimentally 
demonstrates performance of the new 7 method with a laboratory 7 
prototype including the sensorless methods. 

I. Introduction 

Direct torque control of induction machines has 
increasingly become an alternative to field orientation 
methods [1,2]. The classical method of direct torque control 
involves the use of a look-up table to select voltage vectors 
based on torque and stator flux magnitude error [3]. There is 
no current regulator, no pulse width modulation (PWM), nor 
reference frame transformations as in field orientation. With 
appropriately high sample rates this leads to fast torque 
response and low ripple. In addition, at operating conditions 
where the stator flux vector can be estimated accurately from 
the terminal voltage and current, this technique is also 
position and speed sensorless. However, at extremely low 
and zero speeds, the sensorless implementation of this 
technique suffers the same performance degradation as any 
control technique based on the estimate of stator flux using 
only fundamental voltage and current. 

An alternative method of direct torque control is based 
on the deadbeat (inverse) solution to the machine equations 
[4-7]. The deadbeat solution is similar to the classical direct 
torque control method in that it controls torque and stator flux 
directly, without an intermediate current loop. It is different, 
however, in the calculation of the voltage vector to be applied 
to the machine. In the deadbeat solution, an inverse model is 
used to calculate the theoretical voltage vector needed to 
move the machine torque and stator flux to the desired values 
in one sample period. This voltage vector is then synthesized 
over the sample period by the use of PWM modulation 
techniques. However, the calculation of the voltage vector 
requires the solution of a quadratic equation with several 
parameter dependent coefficients. Insight into the operation 
of the control is lost with a purely algebraic approach to the 
solution of the quadratic equations. 


This paper presents a new direct torque control strategy 
where stator and rotor flux are chosen as state variables in the 
deadbeat solution [8,9]. The use of stator and rotor flux as 
state variables, represented in the stator flux oriented 
synchronous reference frame, allows the construction of an 
intuitive graphical depiction of the necessary voltage vector 
to achieve the commanded torque and stator flux magnitude 
values in one time step. The graphical depiction changes as 
operating conditions or parameters vary; thus the resulting 
change in the necessary voltage vector can clearly be seen. 
Conversely, the impact of the selection of the wrong voltage 
vector on both the stator flux and torque errors can also be 
seen. In addition, practical operating limits can be shown on 
the same graph, thus presenting a good visualization of the 
current and voltage limitations. It is further shown that the 
structure of this control strategy is suitable for use with the 
self-sensing position estimation technique [10]. This enables 
low 7 and zero speed position and speed sensorless direct 
torque control. 


II. Formulation of a Stator-Rotor Flux 
Deadbeat Control Algorithm 


The state equations for the induction machine, using the 
stator and rotor flux as state variables in the stationary 
reference frame, are as follows [11]. 


P^qds “ Vqds ' 


♦(; 


R^Lm^ - 

OLgLy 


l^qdr 


P^qdr -^- CTLr +j®tj^qdr + ( Lr cL s )^ ds 

3PLm 

= 4cL s Lj. ^ ds x ^qdrl 


(i) 


( 2 ) 


(3) 


A discrete time form of (l)-(3) is shown in (4)-(6) that is 
valid for small values of the sample time, t s , for which the 
rotor speed, , changes negligibly [8,9]. 

^qds(^+ 1 ) - ^qdsOO = 

v qds (k)t s - + | a L s L r ^ dr<li),s (4) 

^qdrO c+ 0 ‘ ^qdr(K) = 

+ jtA-^^qdrf^)Is + ^7 ^ ^qds(k)t s (5) 

3PLm 

T e (k+ 1 ) = { ^qdsfk+ 1 ) X Xqdr(k+ 1 ) } (6) 
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Equations (4)-(6) can be combined to form and expression 
for the change in torque, AT e (k) = T e (k+1) - T e (k). If the 
d-axis of the excitation reference frame is aligned with the 
stator flux and the terms proportional to tf are neglected, a 
very useful relationship results as shown in (7) [8,9]. 


ATe(k) _ M 3PL m 

h ~ Te< k '[aL r ' oL s J + 4aL s 4 


{ -«v A. d e s (k) X£(k) + Vq S (k) X|(k) - V ds (k) Xq r (k) } (7) 


Equation (7) can be rearranged as follows to show the linear 
relationship between Vq S (k)t s and V ds (k)t s values which can 
be used to provide a given value of AT e (k). 


Vqs(k) t s = 


4o'L s L r 

3PL m >4(k) 


AX e (k) + X e (k) 


aL r 



+ (OfX^stk) + V ds (k) t s 


Aqr(k) 

\dr(k) 


( 8 ) 


If (S) is plotted in a d-q plane with V ds (k)t s and Vq S (k)t s as 
the d- and q-axis variables respectively, the voltage loci for a 
given AT e (k) is a straight line. This line is parallel to the 
rotor flux vector, ta^k), as shown in Fig. 1 [8,9]. (Figs, 1-6 
use the convention tq d =fq-jf d .) 



Fig. L Constant AT e (k) line 

One representation of the multiple possible voltage (volt- 
sec) vectors [Vq S (k)t s V ds (k)t s ] that could yield the desired 
change in motor torque, AX e (k), is shown in Fig. 2. 



Fig. 2. Graphical representation of multiple voltage vectors 
for machine torque change AT e 
From (4) the discrete time expression for stator flux, 
neglecting stator resistance, is 

Xq ds (k+1) - Aq ds (k) = v qds (k)t s . (9) 


Equation (9) can also be shown graphically. Fig. 3 shows 
the plot of lXq ds (k+l)l and X ds (k), where X ds (k) = lXq ds (k)!. 
There are multiple voltage vectors, scaled by t s , which will 
move the flux magnitude from a value of X ds (k) to 
l^qds(k+ Ol- 



Fig. 3. Graphical representation of multiple voltage vectors 
for flux magnitude increase 

The flux circle shown in Fig. 3 can be redrawn to show the 
change in flux magnitude, AlXq d s(k)l = lXq ds (h+Of - ^ds(k)> 
where X ds (k) = lXqd s (k)L This circle is centered on the present 
value of stator flux, X ds (k), with a radius equal to the 
magnitude lAq ds (k+l)l. This adjustment means that the 
possible voltage (volt-sec) vector loci shown in Figs. 2 and 3 
now all begin at the origin of the plot. Thus, for a given set 
of AT e (k) and AlXq ds (k)l requirements, the voltage vector 
which will solve both conditions simultaneously can be 
determined from the intersection of the AT e (k) line and the 
Aa qds (k)l circle. This is shown in Fig. 4. 



Fig. 4. Graphical representation of a voltage vector which 
satisfies both AT e (k) and lXq ds (k+l)f requirements 

The solution, V qds (k)t s , to the intersection of the AT e (k) 
line and the AiA,q ds (k)l circle is bounded by the available dc 
bus voltage, the sample time and the inverter current limit. 
The range of voltages that can be synthesized from a two- 
level inverter can be represented as a hexagon in the d-q 
plane [11], Thus the bound for Vq ds (k)t s is a hexagon with 
sides equal to 2/3 V dc t s as shown in Fig. 5. The hexagon is 
shown at a static position in Fig. 5 but it actually rotates at the 
synchronous speed because the figure is in the synchronous 
reference frame. 
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Fig. 5. Graphical representation of voltage and current limits pig ^ Graphical representation of operating regions for 

on the possible voltage vector solutions steady state conditions 


Fig. 5 also shows the maximum and minimum values of 
AT e based on a steady state solution using the inverter current 
limit, liqdsmax^ anc * the present value of T e (k). The maximum 
torque can be calculated as a function of the stator flux 
magnitude and the maximum stator current as follows. For a 
stator flux oriented system, the torque can be expressed as 

T e = X X d S iqs ( 10 > 

In terms of current magnitude, this becomes 



(ID 


The variables i ds and can be related through the slip 
frequency by using the following two relationships which are 
given in [11]. 


^qds “ Laslqds + ^qdr 

( ~Rr . . ^ _ /L m R r ^ t 

P^qdr = L ^qdr) + I j l qds 


( 12 ) 

(13) 


Assuming stator flux orientation and steady state conditions, 
(12) and (13) can be used to form an expression for i ds in 
terms of the stator current magnitude and the stator flux. 


Ms ' 


LqsWiqds ^ + ^ds 

^Ml+a) 


(14) 


Equation (14) can be substituted into (11) to find a condition 
on the maximum torque that is possible for a specified 
maximum stator current magnitude and stator flux value. 



LgsTs^qdsmax^ + 

^ds^sO - ^) 


(15) 


Finally, using (15) and the present value of torque, the 
maximum and minimum change in torque can be calculated. 


AT emax = +T em ax - T e (k) (lb) 

ATemin = "Temax “ T e (k) (12) 


It is interesting to consider the plot for a steady state 
condition (AX e (k)=0, AiA^sllOMO). This is shown in Fig. 6. 


Any AX e lines to the right of the AT e (k)=0 line represent an 
increase in torque. Any AT e lines to the left of AT e (k)=0 
represent a decrease in torque. Any AtXq ds (lc)l circles within 
the AlXq ds (k)l=0 circle represent a decrease in stator flux 
magnitude. Any AlXq ds (k)l circles outside Al^ ds (k)i=0 
represent an increase in flux. Thus it is seen that there are 
four operating regions within the hexagon as shown: 

1 . Increase torque, decrease flux 

2. Increase torque, increase flux 

3. Decrease torque, decrease flux 

4. Decrease torque, increase flux 

The standard table look up method of direct torque control 
allows only the voltage vectors represented by the vertices of 
the hexagon (plus the two zero vectors) to be selected for the 
duration of the sample period. However, if a. pulse width 
modulation technique is used, an average voltage vector over 
the sample period can be synthesized which lies anywhere 
within the hexagon. 

A control algorithm can be developed based on (8) and (9) 
and using Fig. 4 if AX e (k) is set equal to X e (k+1)* - T e (k) and 
AlXqdsOOl = iXq ds (k+l)l* - >td e s(k). The voltage vector 
calculated by the control algorithm is the intersection of the 
AT e (k)* line with the A0i qds (T)l* flux magnitude circle which 
is the solution of (8) and (9). The voltage vector can be 
synthesized using space vector modulation techniques to 
calculate the inverter switch duty cycles [12]. Fig. 7 shows 
the block diagram of the control system in which the 
proposed algorithm was implemented. (The additional 
voltage command, Vhf ss (k), is only necessary for the self- 
sensing position and speed estimation as explained in Section 
V.) 

The control sequence is as follows. 

1. Is AT e (k)* within the limits of At emax and AT em j n ? W 
not, AT e (k)* is set equal to the closest limit. 

2. Do the AT e (k)* line and the AiX qds (lc)r circle intersect? 
If so, the voltage vector, V q ds(k) is calculated. If not, 
T e (k+1)* is reduced until the corresponding AT e (k)* 
line is tangential to the ADi qds (T)l* circle. Then the 
voltage vector to this point is calculated. 
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Fig. 7. Proposed deadbeat control block diagram 

3. Is the desired voltage vector within the hexagon? If 
so, space vector modulation is used to calculate the 
corresponding inverter switch duty cycles. If not, the 
magnitude of V qds (k) is reduced until it lies on the 
hexagon boundary. Then space vector modulation is 
used to calculate the inverter switch duty cycles. 
Alternate ovemiodulation strategies could also be 
used — for example the most appropriate hexagon 
comer could be selected based on the desired change 
in torque and flux which would be equivalent to using 
standard direct torque control. 

in. S imulation Results 

The experimental motor used in this paper is a specially 
designed high speed (rated speed just over 23,000 rpm) 
induction machine for a NASA electro-mechanical actuator 
research project. The parameter values are given in the 
Appendix. The simulation conditions are set to match the 
experimental conditions. Experimentally, the machine was 
limited to 10,000 rpm or less due to concern about the 
condition of the bearings. All of the simulations neglect the 
PWM switching harmonics. 

The controller was analyzed in two ways. First, to 
investigate its small signal stability, the operating point 
model was formed using the Matlab™ LTI function. To form 
the operating point model, the mechanical dynamics were 
neglected (constant speed) and perfect flux estimation was 
assumed (the flux observer was not included in the operating 
point model). Pole-zero migration plots were then created for 
a range of speeds from low^ speed to rated speed (180 rpm to 
23,000 rpm) at the rated torque and rated flux operating point. 

The second analysis was based on a time domain, nonlinear 
model simulation of the complete system including the flux 
observer and the mechanical dynamics. The speed, torque 
and stator flux responses to a square wave torque command 
are shown. 


Fig. 8 shows the pole-zero migration plot from the LTI 
analysis for a 100 jttsec sample time (used experimentally). 
The result is seen to approximate the expected deadbeat 
response but there is not exact pole-zero cancellation and the 
free pole is not exactly at the origin. This improves as a 
smaller and smaller sample time is used. (The NASA motor 
has relatively small time constants as can be seen in the 
Appendix.) 



Fig. 8. Pole-Zero plot for 100 psec sample time and no 
calculation delay 


The time domain results shown in Figs. 9 and 10 
demonstrate good performance. The torque response is seen 
to be a square wave as expected. Fig. 10 also shows the 
effect of neglecting the resistance in deriving (9). For a stator 
resistance equal to zero, the stator flux is at the commanded 
value of 1 per unit. However, the stator flux magnitude is 
slightly reduced from the commanded value for a stator 
resistance of , 1 8 ohms. 




Fig. 9. Torque and speed 
response for a 100 psec 
sample time and no 
calculation delay 


Fig. 10. Stator flux 
magnitude for 100 psec 
sample time, no calculation 
delay and two values of 
stator resistance 


In the ideal case, it Is assumed that the command values For 
the manipulated variables (V qs *(k), V ds *(k)) are calculated in 
zero time based on the present sampled value of the measured 
feedback variables (i qs (k), i ds (k)). Realistically, however, 
there will be a calculation delay between when the feedback 
variables are measured and when the commanded values for 
the manipulated variables are updated because the required 
microprocessor calculations can not be done in zero time. 
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The command values for the manipulated variables Vqs(k) 
and V d *(k) ultimately become pulse width commands for the 
gate drives in the inverter. In this implementation, the 
sampling time is synchronized with the PWM generation. 
Thus the PWM timers are updated at the same time as the 
feedback variables are measured. This results in a one step 
time delay for the control because the PWM timers are 
updated based on the calculations from the feedback data of 
the previous sample. 

The effect of this one step time delay can be seen in the 
pole-zero plot of Fig. 1 1 . The free pole and the pole-zero 
cancellation pair that were at or near the origin in Fig. 8 have 
now moved to the edge of the unit circle. The time domain 
simulations also show a more oscillatory response as seen in 
Figs. 12 and 13. 



Fig. 11. Pole-Zero plot for 100 jisec sample time and one 
step time delay 



Fig. 12. Torque and speed 
response for a 100 psec 
sample time and a 1 time 
step calculation delay 



Figure 13: Stator flux 
magnitude for 1 00 psec 
sample time and a 1 time 
step calculation delay 


For a similar deadbeat type of controller, Ref [5] shows 
that a modification to AT e (k)* and AlXq ds (k)r will move the 
system poles to lie further within the unit circle and reduce 
the oscillations. A similar derivation can be shown here. To 
include the one-step time delay, (7) can be expressed as 
follows. 


T e (k+2)-T e (k+l) 


aLr 


R s 

= Te(k+l)|-r^-— 1 + 


3PL m 
I 4ctL s L t 


{-^(k+D^k+l) + Vqs(k)* ^ d e r (k+l) - v d e s (k)* A£(k+i)} 

(18) 


From (8) an expression for the control variables is Vq S (k)* 
and V ds (k)* is 


% ^ ( k ) 

V q Wts-V d W‘ s 7f— = o^dsflOts + 




4 <7 1 - S L; 
3PL„Aj e r (k) 


t e (k + ir-x e (k ) + t e (k)i— + — It s n9) 


Rr 


Under constant stator flux operation and assuming that 
^qdi-ik+1) ~ Xqdrfk), (18) and (19) can be combined to form 
an approximate relationship between the commanded torque 
and the actual torque as shown in (20). 


T e (k+2) - T e ( k+ 1 ) = t e (k+ 1 ) 


+ T e (k+ 1 )* 


( ~R r _M 
[ah, ' oL s J ts 

f Rr ^ R s 

I aLj- aL s 


- Te(k) + T e (k) 


( 20 ) 


From (20), the characteristic equation of this simplified 
transfer function is 


{'■(£*£)•>('■(£*£) 


ts =0. (21) 


For small values of t s , these poles can be seen to lie close 
to the unit circle boundary. Ref [5] suggests modifying the 
control equations as shown in (22) and (23) (0<C<1) to 
reduce the oscillations and potential instability problems. 


AT e (k)* = C (T e (k+1)* - T e (k)) 


( 22 ) 


AlXqds(k)l* = c 0^ds(k+l)r - l^dsOOO (23) 


Using (18), (19) and (22), the characteristic equation of the 
modified transfer function is 



(24) 


Ref [5] studies the effect of various values of C on the 
response of the system. In the present w ork, it was found that 
C=.8 yielded a fast response with no instability. The 
corresponding LTI pole-zero migration plot and time domain 
responses are shown in Figs 14-16. 



Fig. 14. Pole-Zero plot for 100 psec sample time, one step 
time delay and C=.8 
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and 20 (without deadtime) to Figs 17 and 18 (with deadtime). 
It is seen that without deadtime compensation, both the 
torque and stator flux magnitudes were reduced. 


I 

i 


0.2 0 4 0.6 0.6 1 

time, wcorKfc 

Fig. 16. Stator flux 
magnitude for a 100 psec 
sample time, I time step 
delay and C= 8 

IV. Experimental Results 

The experimental test set up consisted of a dc power 
supply, an intelligent power module, a dSpace 1 103™ digital 
controller and the NASA test motor. The results were 
captured using dSpace software. The data files were then 
plotted using Matlab . Twx> phase currents were measured 
and an encoder was used for position feedback. The position 
feedback was necessary because a flux observer based on the 
current model was used. (In Section V, the position feedback 
is provided by the self-sensing algorithm and the encoder 
information was used for comparison purposes only.) The 
speed was calculated in the controller using the position 
information. The torque and flux were estimated in the 
controller using current and position information. 

Figs. 17 and 18 show the torque, stator flux and speed for 
the same conditions as in Figs. 15 and 16 in simulation. 



Fig. 15. Torque and speed 
response for a 100 psec 
sample time, 1 time step 
delay and C=.8 




Fig. 17. Experimental results: 
torque and speed response 
for a 100 psec sample time, 
one time step delay and C=.8 


Fig. 18: Experimental 
results: stator flux magnitude 
for a 100 psec sample time, 
one time step delay and C=.8 


In general, the response is as predicted. In the actual 
implementation, there is always a deadtime, or blanking time, 
in the inverter so that the two switches across a leg do not 
conduct at the same time. In addition, there is a voltage drop 
across the conducting devices. Both of these effects result in 
a lower voltage being applied to the machine than is actually 
commanded. This "deadtime voltage" can be calculated [13], 

Vdeadrime^ (25) 

Vd is the magnitude of the voltage due to the combination 
of the deadtime losses and the conduction losses and k is the 
sector of the d-q plane in which the current vector is located, 
k=l,2...6. This controller was found to be sensitive to 
deadtime compensation as can be seen by comparing Figs. 19 



Fig. 19. Experimental results: 
torque and speed response 
for a 100 psec sample time, 
one time step delay and 
t s '=l 25 psec, no deadtime 
compensation 


Fig. 20, Experimental results: 
stator flux magnitude for a 
100 psec sample time, one 
time step delay and t s '=125 
psec, no deadtime 
compensation 


V. Low Speed Sensorless Implementation 

To date, the self-sensing method has been demonstrated 
only with the field orientation method of control. The use of 
self-sensing in a deadbeat direct torque control offers the 
potential of full speed range sensorless operation with the 
flux estimate based on the back emf method at higher speeds 
and on the self-sensing position estimate at lower speeds. 

To estimate the rotor position angle, 0 r , the self-sensing 
technique requires a machine with a magnetic saliency related 
to the rotor position [10], The NASA motor rotor was 
designed to produce a position dependent magnetic saliency 
by changing the shape of the rotor bars as a function of 
position as shown in Fig. 21. (The rotor is 1.3 inches in 
diameter). 



Fig. 21. Rotor cross-section of NAS A machine 

Additionally, the self-sensing technique requires a high 
frequency carrier signal. The structure of a deadbeat type of 
direct torque control lends itself to the use of the self-sensing 
method of position estimation because the necessary high 
frequency signal can be easily added to the fundamental 
voltage command as showm in Fig. 7. Thus the com man ded 
voltage into the pwm modulator will consist of the 
fundamental voltage command, V^Ck) and the high 
frequency voltage, Vjjf ss (k). The resulting current is filtered 
as described in [14] and a signal at twice the rotor position is 
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generated due to the position dependent magnetic saliency. 
This signal can be tracked in a closed loop, saliency image- 
tracking observer to produce position, velocity, acceleration, 
and disturbance torque estimates. 

Figs. 22 and 23 show the results of a no load closed loop 
low speed sensorless control using the self-sensing position 
estimate as feedback for the controller. Fig. 22 shows the 
spectra of the negative sequence current (as defined in [14]) 
for a constant 3 Hz speed. The component due to the rotor 
saliency is clearly visible at 6 Hz. Fig. 23 shows a no load 
speed reversal from -3 Hz to + 3 Hz. 



Fig. 22. Negative sequence Fig. 23. Sensorless speed 
current spectra for constant reversal of proposed 

speed 3 Hz operation controller with estimated 

rotor position and speed 
feedback based on self- 
sensing method 


VI. Conclusions 

This paper has developed a new, deadbeat direct torque 
control method based on stator and rotor flux as state 
variables. This choice of state variables allows a clear 
graphical visualization of the voltage vector solution and the 
inverter operating limits. 

The implementation of the proposed controller was 
evaluated experimentally and found to produce good results. 
The implementation issues which could limit performance 
were also evaluated. 

The controller was found to be sensitive to the one step 
time delay in the experimental implementation. The 
deadtime voltage drop, without appropriate compensation, 
was also found to reduce the torque and flux in the machine. 

The structure of the proposed controller allows the addition 
of a high frequency voltage vector to the commanded 
fundamental voltage vector. This allows the self-sensing 
method of position and speed estimation to be used thus 
demonstrating low, including zero speed, sensorless control. 
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Appendix 

NASA Motor: 2 pole induction, 96 V LN . 400 Hz, 

1^=1. 9e-3 H., L ls =1.25e-4 H., L ir =1.25e-4 H., R r =.105 Cl., 

R s =.09 £L, J=1.02e-4 kg-m^, T en =l N-m, 0^23,030 rpm 
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